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ABSTRACT 

The  joint  detection  of  all  users  in  a  multiple  access 
communication  system  is  known  to  greatly  enhance  the 
system  performance.  The  optimal  joint  detector  has  a 
complexity  which  is  exponential  in  the  number  of  users 
and  is,  therefore,  impractical.  Suboptimal  algorithms 
attempt  to  achieve  near  optimal  performance  with  re¬ 
duced  computational  complexity.  The  suboptimal  joint 
detection  literature  deals  only  with  receiver  design,  al¬ 
lowing  for  the  use  of  any  set  of  signature  waveforms. 
Using  the  structure  of  the  wavelet  packet  transform  [1] 
we  are  able  to  choose  signal  sets  and  design  a  detection 
algorithm  which  results  in  significantly  lower  computa¬ 
tional  complexity  than  other  proposed  suboptimal  joint 
detectors.  The  wedding  of  detector  development  and 
wavelet  packet  signal  set  design,  therefore,  becomes  an 
interesting  problem.  Convergence,  performance  and 
receiver  complexity  of  suboptimal  MA  detectors  cru¬ 
cially  depend  upon  the  signal  set  structure;  it  is,  thus, 
of  paramount  importance  to  study  the  convergence  of 
our  iterative  algorithm  for  our  signal  structure.  In  this 
paper,  sufficient  conditions  for  the  convergence  of  an 
iterative  joint  detection  scheme  are  developed. 

1.  INTRODUCTION 

A  multiple  access  communication  system  will  typically 
support  a  large  number  of  users  over  a  given  channel. 
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By  allowing  many  users  to  transmit  simultaneously  in 
the  same  frequency  band,  we  are  inducing  interference 
among  users.  This  inter-user  or  multiple  access  inter¬ 
ference  (MAI)  will  degrade  the  performance  of  conven¬ 
tional  receivers,  and  is,  in  fact,  the  major  limiting  fac¬ 
tor  when  the  number  of  users  exceeds  a  threshold.  In 
the  current  literature,  much  attention  is  being  paid  to 
overcoming  the  MAI  effects,  but,  the  signal  design  as¬ 
pect  and  the  receiver  design  aspect  of  the  MA  system 
are  always  considered  separately. 

Assuming  that  a  single  bit  is  transmitted  by  each 
user  within  a  single  time  window,  the  received  signal, 
r,  which  is  the  aggregate  of  each  user’s  transmission 
(plus  noise),  can  be  written  as 

K 

v  =  '^2  bkWk  +  n  =  Wb  +  n,  (1) 

k= 1 

where  K  is  the  total  number  of  users,  bk  is  either  a  1 
or  —1  corresponding  to  a  1  or  0  information  symbol 
of  the  kth  user,  w k  is  the  kth  user’s  discrete  waveform 
represented  as  a  vector,  n  is  an  i.i.d.,  zero  mean,  Gaus¬ 
sian  noise  vector,  and  the  following  definitions  are  used: 
b  =  [&!  &2  '  • '  hx]T,  W  =  [wi  w2  ■  •  •  wjf].  We  assume 
that  all  relative  time  shifts  are  known  and  incorporated 
in  Wfc. 

The  conventional  approach  to  receiver  design  is  to 
use  a  matched  filter, 

h  =  sgn[  w|r],  (2) 

where  sgn  represents  the  signum  functional  and  bk  de¬ 
notes  the  estimated/detected  symbol  of  the  kth  user. 
The  conventional  matched  filter  (2)  represents  the  op¬ 
timal  receiver  in  the  case  where  the  MAI  is  assumed 
normal  and  white.  Under  this  assumption,  as  users  are 
added  to  the  system  the  MAI  raises  the  noise  level;  this, 
in  turn,  limits  the  matched  filter  performance  and,  ul¬ 
timately,  the  reliable  throughput.  For  this  reason,  the 


signal  design  aspect  of  MA  is  typically  restricted  to  a 
general  attempt  at  minimizing  correlation  among  user 
waveforms. 

The  MAI,  however,  is  not  an  additive  white  Gaus¬ 
sian  noise  process;  it  possesses  a  great  degree  of  struc¬ 
ture  which  can  be  exploited  in  achieving  a  greater  trans¬ 
mission  capacity  (e.g.  more  simultaneous  users)  than 
the  system  which  uses  the  conventional  detector.  The 
optimal  joint  detector  (which  accounts  for  the  MAI) 
maximizes  the  log-likelihood  function  [2]  and  results  in 

b  =  arg  [  max  2rrWb  -  (Wb)TWb  ].  (3) 

b:  {i,  -i}“ 

This  maximization  is  over  2K  possible  b  vectors  (a  com¬ 
putational  complexity  that  is  exponential  in  I(,  the 
number  of  users).  This  optimal  method  offers  sub¬ 
stantially  higher  performance  than  the  conventional  de¬ 
modulator  (2)  (e.g.  it  will  afford  a  relatively  larger 
number  of  users  successful  utilization  of  an  MA  sys¬ 
tem),  but  is  of  little  practical  use  on  account  of  its 
computational  burden  [3]. 

Recent  communication  literature  addresses  the  gen¬ 
eral  notion  of  suboptimal  joint  detection  which  offers, 
with  minimal  performance  loss,  computational  improve¬ 
ment  relative  to  the  optimal  method  while  achieving  a 
significant  improvement  over  the  conventional  detec¬ 
tor  [4,  5,  6j.  These  approaches  only  address  the  re¬ 
ceiver  design  problem.  The  convergence,  performance, 
and  receiver  complexity  of  such  suboptimal  MA  detec¬ 
tors,  however,  crucially  depend  upon  the  structure  of 
the  signal  set  W.  In  light  of  this,  our  work  in  MA 
systems  couples  detection  and  waveform  design. 

Our  method  for  waveform  design  takes  full  advan¬ 
tage  of  the  hierarchical  nature  of  the  wavelet  packet 
transform  to  control  MAI  by  providing  structure  that 
can  be  exploited  in  receiver  design.  Our  waveforms 
possess  a  nested-type  of  correlation  structure  which  is 
well  suited  for  iterative  bit  estimation.  As  shown  in  [7] 
by  using  our  wavelet-packet-based  signature  waveforms 
in  lieu  of  the  well  known  direct  sequence  spread  spec¬ 
trum  waveforms,  a  great  reduction  in  computational 
complexity  of  other  proposed  detectors  results. 

2.  THE  JOINT  DETECTOR 

In  the  MA  scenario  described  in  Equation  (1),  we  ob¬ 
serve  the  signal 

r  =  Wb  +  n.  (4) 

We  wish  to  derive  an  iterative  algorithm  to  estimate 
b.  We  know  that  b  6  { 1 ,  —  1 } /v  and  we  use  as  the  first 
estimate  the  result  from  the  matched  filter  detector  (2). 
Denoting  the  first  estimate  of  the  bit  vector  as  b(l),  we 
have 

b(l)  =  s^n[WTr],  (5) 


Without  noise,  i.e.  n  =  0,  we  have  r  =  Wb  and  our 
first  estimate  would  be 

b(l)  =  s<i»[WTWb].  (6) 

For  obvious  reasons  we  require  the  true  bit  vector,  b,  to 
be  a  fixed  point  of  our  iteration,  at  least  in  the  noiseless 
case.  In  general,  b  ^  WrWb.  For  the  special  case 
in  which  b  =  WTWb,  the  user  waveforms  would  be 
orthogonal  and  the  matched  filter  would  suffice.  In  the 
general  case  (in  which  users  are  correlated)  we  must 
add  (b  -  WTWb)  to  the  argument  of  Equation  (6)  so 
that,  in  the  noiseless  case,  we  have 

b  =  sgn[ WrWb  +  b  -  WTWb] .  (7) 

This  leads  to  our  iterative  detector 

b(m  +  1)  =  sgn[ WTr  -I-  b(m)  —  WrWb(m)],  (8) 
or 

bk(m  +  1)  =  sgn[vTwk  -  ]Pw*w jbj(m)],  (9) 

where  m  denotes  the  iteration  number  and  b(0)  =  0 
so  that  Equation  (8)  reduces  to  the  matched  filter  es¬ 
timate  (5)  for  b(l). 

We  offer  a  simple  interpretation  of  our  joint  detector 
by  looking  at  Equation  (8)  for  noiseless  reception, 
r  =  Wb, 

b(m  +  1)  =  sgn[b(m)  +  WTW(b  —  b(m))].  (10) 

At  each  iteration,  we  update  the  previous  estimate  with 
a  correction  term  which  is  based  on  the  difference  be¬ 
tween  the  actual  bit  vector  and  the  previous  bit  esti¬ 
mate. 

The  joint  detection  procedure  defined  by  Equation  (9) 
is  similar  to  a  more  general  algorithm  proposed  by 
Varanasi  and  Aazhang  [4]. 

Commonly,  MA  systems  employ  direct  sequence  spread 
spectrum  (DSSS)  or  pseudo-noise  spread  spectrum  (PNSS) 
user  waveforms  (sequences)  where  any  single  user  expe¬ 
riences  interference  from  all  other  users  [8].  Figure  1-a 
shows  the  cross  correlation  matrix,  WrW,  for  a  set  of 
68  pseudo-noise  (PN)  waveforms  in  a  64  dimensional 
waveforms  space  (i.e.  64  samples  per  waveform  dura¬ 
tion).  Each  waveform,  w*,,  was  generated  by  a  ran¬ 
dom  number  generator  such  that  wk  €  {—1,  l}64.  We 
make  two  observations:  1)  when  is  a  pseudo-noise 
spread  spectrum  (PNSS)  waveform,  w^w j  is,  in  gen¬ 
eral,  nonzero  for  all  values  of  j  and  k,  therefore,  Equa¬ 
tion  (9)  must  decode  the  bits  of  all  K  users  to  detect 
a  single  user’s  bit;  2)  each  user  experiences  the  same 


degree  of  MAI.  We  note  that,  typically,  the  number  of 
users  in  a  DSSS  system  does  not  exceed  the  number  of 
dimensions  in  the  waveform  space.  Additionally,  Gold 
or  Kasami  sequences  [8]  which  have  been  designed  to 
have  low  cross  correlations  are  typically  used.  We  in¬ 
clude  this  example  of  the  PN  waveform  set  in  order 
to  illustrate  the  general  notion  of  DSSS  and  the  cross 
correlation  structure  of  PN  sequences. 

Figure  1-b  shows  the  cross  correlation  matrix  fox- 
one  wavelet  packet  signal  set  in  which  68  user  wave¬ 
forms  exist  in  a  64  dimensional  waveform  space.  (See  [7] 
for  the  wavelet  packet  signal  set  development.)  In  con¬ 
trast  to  the  PNSS  observations  made  above:  1)  each  of 
our  wavelet  packet  waveforms,  w is  correlated  with 
only  a  fraction  of  the  signal  set;  2)  each  user  experi¬ 
ences  a  different  degree  of  MAI.  We  immediately  see 
that  our  waveforms  offer  two  advantages  in  a  joint  de¬ 
tection  MA  system:  1)  due  to  the  reduction  in  the 
number  of  correlated  users,  a  large  degree  of  computa¬ 
tional  simplification  of  the  joint  detector  is  realized;  2) 
the  hierarchical  structure  of  MAI  experienced  by  each 
user  lends  itself  to  the  iterative  nature  of  our  algorithm. 

In  the  next  section  we  study  the  convergence  of  an 
iterative  joint  detection  algorithm.  We  set  a  framework 
within  which  wavelet  packet  multiple  access  (WPMA) 
signal  optimization  as  well  as  joint  detection  can  be 
considered  coincidentally. 

3.  STUDY  OF  CONVERGENCE 

We  are  interested  in  the  convergence  of  Equation  (10), 
the  recursive  joint  detection  algorithm  for  the  noiseless 
case.  We  rearrange  Equation  (10)  and  incorporate  the 
scalar,  p,  to  allow  for  speedy  convergence. 

b(m  +  1)  =  sgn[b  +  (gWTW  —  I)(b  —  b(m))j.  (11) 

Proposition:  The  optimal  convergence  speed  of  the 
recursion  in  Equation  (11)  is  achieved  for  p  =  ^|g, 
where  A  and  B  are,  respectively,  the  minimum  and 
maximum  eigenvalues  of  the  matrix  WTW.  Further¬ 
more,  the  convergence  is  guaranteed  if 

B  <  3 A. 


Proof:  Let  e(m)  denote  the  error  vector  with  its  jth 
element  defined  by 


ej  (m) 


A 


0  ,  if  bj(m)  =  bj 

1  ,  otherwise. 


The  error  vector  indicates  which  bits  are  incorrect  in 
the  estimate  b(m),  for  m  >  1.  We  set  e(0)  =  0.5 77, 


(a) 


(b) 

Figure  1:  WTW  displayed  as  a  gray  scale  image  where. 
White  corresponds  to  zero  and  black  corresponds  to 
unit  magnitude,  a)  PN  sequences  b)  WPMA  wave¬ 
forms 

where  7?  =  [1, . . . ,  1]T.  Letting  B  =  diag{b}  allows  us 
to  express  the  error  in  our  estimate  as 

(b  —  b(m))  =  2  Be{m). 

Let  Z  =  (pWTW  —  I)  so  that  we  may  rewrite  equation 
(11)  as 

b(m  +  1)  =  sgn[£>?7  +  2ZSe(m)] 

=  B  sgn[r]  +  2BZBe(m)],  (12) 

where  we  used  the  equality  BB  =  I.  From  Equa¬ 
tion  (12),  we  see  that  an  error  is  made,  i.e.  bk(m  +  1) 
is  incorrect  and  efc(m  +  1)  =  1,  if 

[2BZBe(m)]k  <  — 1. 

In  more  general  terms,  if 

||  2BZBe(m)  ||2<  p  +  1,  p  €  {0, 1, 2,  •  •  •} 


||  2BZBe(m)  ||<  \/p+  1, 

then  b (m  +  1)  has  at  most  p  +  1  incorrect  bits  and 
||  e(m  +  1)  ||<  \Jp  +  1. 

Now,  if  b(m  +  1)  has  exactly  p  +  1  incorrect  bits, 
then 

||  e(m  +  1)  ||=  y/p+l 

and  we  have 


||  2BZBe(m  +  1)  ||  <  2  ||  BZB  ||  ||  e(m  +  1)  || 

=  2  ||  BZB  ||  VpT T.  (13) 

In  the  next  iteration  of  the  estimate,  b(m  +  2),  we 
tolerate,  at  most,  p  incorrect  bits,  i.e. 

||  e(m  +  2)  ||<  y/p+1, 

so  that  the  number  of  bit  errors  decreases  at  each  step 
of  the  recursion 


...  <||  e(m  +  2)  ||<||  e(m  +  1)  ||<  . . . 

and  convergence  is  insured.  Hence,  for  convergence,  we 
need 

II  BZB  ||  <  i. 

In  particular,  in  order  to  have  optimal  convergence 
speed,  we  must  minimize  the  norm  of  (BZB).  We  refer 
here  to  the  matrix  norm 

II  Cr  II 

II  C  ||=  max  V-/  =  \X(C)\maXi 

where  A (C)  denotes  the  eigenvalue  of  the  matrix  C.  It 
is  easy  to  show  that  \(BZB)  =  A(Z)  and 


||  Z  ||=  \X(Z)\max  =  max(|pA  -  1|,  \pB  -  1|), 

where  A  =  Amm(WTW)  and  B  =  Amoa;(WTW), 

0  <  A  <  B  <  oo.  Additionally,  the  norm  of  Z  is 
minimal  if  we  choose  p  =  so  that 


It  follows  that  the  recursion  in  (12)  is  guaranteed  to 
converge  if 

11  z  ll=  <  1/2  ^  S<3A  n 

The  convergence  conditions  derived  above  are  suf¬ 
ficient,  only.  In  [7],  we  arbitrarily  chose  WPMA  signal 
sets  for  which  A  =  0,  i.e.  the  user  signature  waveforms 


comprise  a  linearly  dependent  set,  but  convergence  is 
observed  in  Monte  Carlo  simulations  using  the  itera¬ 
tive  detector  of  Equation  (10).  The  goal  of  future  work 
is  to  find  both  necessary  and  sufficient  conditions  on 
W  to  guarantee  convergence  of  a  joint  detection  al¬ 
gorithm.  These  conditions  will  elucidate  modifications 
of  the  joint  detector  so  that  optimal  performance  is 
achieved  for  any  given  number  of  users  communicating 
within  a  given  time-bandwidth  space. 

4.  CONCLUSION 

The  importance  of  joint  detection  of  all  correlated  users 
in  a  multiple  access  communication  system  motivates 
our  research  in  the  area  of  the  coincidental  design  of 
user  waveforms  and  joint  detection  algorithms.  The  hi¬ 
erarchical  correlation  structure  of  wavelet-packet-derived 
signals  appears  to  be  well  suited  to  the  development  of 
iterative  joint  detectors. 

In  this  paper,  we  set  a  framework  within  which 
wavelet  packet  multiple  access  signal  optimization  as 
well  as  joint  detection  can  be  considered  together.  We 
derived  one  such  iterative  joint  detector  and  developed 
sufficient  conditions  for  its  convergence. 
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